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ABSTRACT 

We consider a colocated MIMO radar scenario, in which the 
receive antennas forward their measurements to a fusion cen- 
ter. Based on the received data, the fusion center formu- 
lates a matrix which is then used for target parameter esti- 
mation. When the receive antennas sample the target returns 
at Nyquist rate, and assuming that there are more receive an- 
tennas than targets, the data matrix at the fusion center is low- 
rank. When each receive antenna sends to the fusion center 
only a small number of samples, along with the sample index, 
the receive data matrix has missing elements, corresponding 
to the samples that were not forwarded. Under certain con- 
ditions, matrix completion techniques can be applied to re- 
cover the full receive data matrix, which can then be used in 
conjunction with array processing techniques, e.g., MUSIC, 
to obtain target information. Numerical results indicate that 
good target recovery can be achieved with occupancy of the 
receive data matrix as low as 50%. 

Index Terms — Array processing, compressed sensing, 
matrix completion, MIMO radar, MUSIC 

1. INTRODUCTION 

Multiple-input and multiple-output (MIMO) radar systems 
have received considerable attention in recent years due 
to their superior target estimation performance. Colocated 
MIMO radar systems exploit waveform diversity to formu- 
late a long virtual array with number of elements equal to the 
product of the number of transmit and receive antennas. As a 
result, they achieve higher resolution than traditional phased 
array radars for the same amount of data [EQlEO. Compressed 
sensing (CS) enables MIMO radar systems to maintain then- 
advantages while significantly reducing the required mea- 
surements per receive antenna 00- In CS-based MIMO 
radar, target parameters are estimated by exploiting the spar- 
sity of targets in the angle, Doppler and range space, referred 
to as the target space. For CS-based sparse target estimation, 
the target space needs to be discretized into a fine grid, based 
on which the CS sensing matrix is constructed. However, per- 
formance of CS-based MIMO radar degrades when targets 
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fall between grid points, a case also known as basis mismatch 
in the CS literature. 

Another approach related in spirit to CS is that of matrix 
completion. Matrix completion aims to recover a low-rank 
data matrix from partial samples of its entries by solving a re- 
laxed nuclear norm optimization problem [6|[7|. Array signal 
processing with matrix completion has been studied in [ 8 1 [ 9 ] . 
To the best of our knowledge, however, matrix completion has 
not been exploited for target estimation in colocated MIMO 
radar. Our paper is related to the ideas in [9 1 in the sense 
that matrix completion is applied to the received data matrix 
formed by an array. However, due to the unique structure of 
the received signal in MIMO radar, the problem formulation 
and treatment in here is different than that in |0. 

The main idea of our work is as follows. We consider 
a colocated MIMO radar scenario in which receive antennas 
forward their measurements to a fusion center. Based on the 
received data, the fusion center formulates a matrix, which is 
then used for estimating the target parameters. When the re- 
ceive antennas sample the target returns at Nyquist rate, and 
assuming that there are more receive antennas than targets, 
the data matrix at the fusion center is low-rank. When each 
receive antenna sends to the fusion center only a small num- 
ber of samples, along with the sample index, the receive data 
matrix has missing elements, corresponding to the samples 
that were not forwarded. Under certain conditions, matrix 
completion techniques can be applied to recover the full re- 
ceive data matrix, which can then be used in conjunction with 
parametric methods such as MUSIC to obtain target informa- 
tion. Compared to CS MIMO radar, our proposed method has 
the same advantage in terms of reduction of samples needed 
for accurate estimation but it avoids the basis mismatch issue 
inherent in CS-based approach. 

The rest of the paper is organized as follows. The colo- 
cated MIMO radar system model is described in Section |2] 
Background of noisy matrix completion is introduced in Sec- 
tion 3. The applicability of matrix completion to MIMO radar 
is discussed in Section 4, and numerical results are given in 
Section 5. Finally, Section 6 provides concluding remarks. 

2. COLOCATED MIMO RADAR SYSTEM 

We consider a MIMO radar system that employs colocated 
transmit and receive antennas, shown in Fig. [TJ We use 
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Fig. 1. Colocated MIMO radar system under the ULA model. 



ter, received data during the qth pulse can be written as ifTUl 

X ? =B(fl)SD,A T (0)S + W, = Z„+W 9 , (3) 
where S = diag Qfii, Pk\)\ D 9 = diag(d g ), with 
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S = 



[s (0T S ) , s ((L — 1) T s )]; and W g is a Gaussian noise 
matrix. Note that both matrices S and D q are rank-X, while 
the rank of matrix S is min {M t , L}. Thus, for M t > K the 
rank of the noise free data matrix Z q = B (0) SD 9 A T (0) S 
is K. In other words, the data matrix Z q is low-rank based on 
the assumption that M r ^> K. 



M t and M r to denote the numbers of transmit antennas and 
receive antennas, respectively. Although our results can be 
extended to an arbitrary antenna configuration, the results 
here are presented for the case in which the transmit and 
receive antennas form a uniform linear array (ULA) with 
inter-element spacing between transmit and receive antennas 
d t and d r , respectively. We further assume d t — d r = A/2, 
where A is the wavelength of the carrier signal. Further, the 
waveforms s k (t) ,k = 1, M t , transmitted from the trans- 
mit antennas are assumed to be narrow-band and orthogonal. 
Now suppose there are K point targets in the far field at 
angles 0fc,fe = 1,...,K, each moving with speed The 
corresponding M t x K dimensional transmit steering matrix 
can be expressed as A (0) = [a (9\) , a (Ok)], where 
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and the M r x K dimensional receive steering matrix B (0) can 
be defined in a similar fashion based on the receive steering 
vectors b(0j,). 

In order to estimate the speed of each target, multiple 
pulses need to be transmitted. Let Q be the number of trans- 
mitted pulses, with Tppj being the pulse repetition interval. 
Assume the target reflection coefficients {(3^} ,k — 1, K 
are complex and remain constant during the Q pulses. For 
slowly moving targets, (2i)T p /\ <C 1, where T p is the pulse 
duration) the Doppler shift within a pulse can be ignored, 
while the Doppler changes from pulse to pulse. 

Under the narrowband assumption the received signal at 
the Ith receive antenna can be approximated as [4 | 
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Heres(t) = [s x (*),... ,s Mt (*)] ■ 

Suppose that each receive antenna samples the received 
signal at rate L/T s and forwards the samples to the fusion 
center. Here, T s is the sampling time, L is the number of 
nonzero samples and we assume L 3> K. At the fusion cen- 



3. MATRIX COMPLETION WITH NOISE 

We now provide a brief overview of the problem of recov- 
ering a rank r matrix M e C" lX ™ 2 based on partial knowl- 
edge of its entries, possibly corrupted by noise, i.e., [Y]^ = 
[M]j - + [E]--, £ ri, where, [E] 4 ,- represents noise and 
£1 is the set of observed entries. This can also be expressed as 
Tn (Y) = Vn (M) + Vn (E), where Vq. represents the sam- 
pling operation. According to [7|, when M is low-rank and 
its singular vectors are sufficiently spread, i.e., both the num- 
bers of zero and large elements in the singular vectors are not 
large, M can be recovered by solving a relaxed nuclear norm 
optimization problem, given by 



s.t. 



||Pn(X- 



Y)|| F <5 



(4) 



where denotes the nuclear norm, i.e., the sum of singular 
values of X, while S > is a constant. 

To test the 'sufficiently spread' requirement, the strong 
incoherence property of matrix M with parameter /i has 
been introduced in (6). Consider the singular value de- 

r 

composition (SVD) of M, i.e., M = VJ /OfcU^v^, and 

fe=l 

define P v = J2 u i u f > p v = E v i v f > and T = 

l<i<r l<i<r 

u i v f ■ When the following conditions are satisfied, the 

l<i<r 

matrix M is said to satisfy the strong incoherence property 
with [i = max (hi, 

Al) For all pairs (a, a') £ [n\] x [m] and (6, b') £ [712] x 
[712], there is /ii > such that 



(5) 
(6) 



where e a is the vector with the ath element equal to 1 and 
others being zero, while l Q=a / indicates that it is equal to 1 
when a = a' and otherwise. 

A2) For all (a, 6) £ [nj] x [712], there exists a constant 
\i 2 > such that |T o6 | < p. 2 ^ 7l2 , where T ab is the (a,b) 
entry of the matrix T. 
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Note that it has been shown in that when the maximum 
(in magnitude) values of the K left and right singular vectors 
are bounded, i.e., 



(7) 



with fis = 0(1), then the strong incoherence property is 
guaranteed with fx < \LB\pr. 

Now define n = max (m, H2) and suppose M satisfies 
the strong incoherence property. Then [6 1 establishes in the 
noiseless case that, by observing N randomly selected entries 
with N > C/j, 2 nrlog 6 n for some constant C, the matrix M 
can be recovered exactly with a probability of at least 1 — n~ 3 . 
Further, [7 1 establishes that, when observations are corrupted 
with white noise [E] ■ that is zero-mean Gaussian with vari- 



ance a , the recovery error is bounded as 

4^/i(2 + p) min (ni,n2)S + 26, where p = 
tion of observed entries. 
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4. MATRIX COMPLETION FOR MIMO RADAR 



Fig. 2. Complementary CDF of mi and m.2. Case I: 
M r = 40, L — 128; (a) left singular vectors, (b) right sin- 
gular vectors; Case II: M r = 1000, L = 128; (c) left singular 
vectors, (d) right singular vectors. 



The left singular vectors of Z q defined in (01 are the eigenvec- 
tors of Z g Zf = HSS ff H H , where H = B (0) SD 9 A T (0). 
The right singular vectors of Z q are the eigenvectors of 
S H H ff HS. Since S is orthogonal, it holds that SS H = I. 
Thus, the left singular vectors are only determined by ma- 
trix H, while the right singular vectors are affected by both 
transmit waveforms and matrix H. 

For the problem considered in this paper, it is difficult to 
determine analytically the behavior of the entries of the left 
and right singular vectors of Z q . Instead, we get an idea of 
the behavior of the maximum values and the parameters that 
affect their spread using simulations. We consider a MIMO 
radar setup in which the target direction of arrival (DOA) an- 
gles are uniformly distributed in [—90°, 90°] and the corre- 
sponding target speeds are uniformly distributed in the range 
[150, 450] m/s. In addition, are following complex Gaus- 
sian distribution and kept unchanged for Q — 10 pulses. The 
pulse repetition interval is Tp m = 1/ 4000 second, and the 
carrier frequency is / = 10 9 Hz, resulting in A = c/f 
I)..'! meter. Two types of orthogonal waveforms are consid- 
ered: Hadamard and Gaussian orthogonal waveforms. Sev- 
eral cases of parameters are verified. Case I: M r = 40, 
L = 128; Case II: M r = 1000, L = 128; Case III: M r = 40, 
L = 1024. Each case runs for 300 iterations. 

Let mi and to 2 denote the maximum element value of K 
left and right singular vectors of Z q . The complementary cu- 
mulative distribution function (CDF) curves of mi and m.2, 
i.e., Pr (m > rrii) ,i = 1, 2 are plotted in Fig. |2]for Cases I 
and II. It can be seen from these plots that as M r increases, the 
bounds of mi for both K = 2 and K = 10 decrease, while 
the distribution of m,2 does not significantly change when L 



is fixed. Figure [2] (c) shows that as M r gets large, mi gets 
bounded by a small number with high probability. Space lim- 
itations prevent us from displaying more figures of Case III, 
which show that 1112 is bounded by a small number with high 
probability as L gets large. 

Extensive simulations also show that the bounds on mi 
and rri2 scale as 1 j \JM r and 1 / \fL, respectively, with some 
constant J]Ib. For Gaussian orthogonal waveform and K — 
2, Hb ~ 2.4 for the bound of mi (see Fig. [3(a)) and /is w 
6.5 for the bound of 777,2 (see Fig. [3] (b)). For K = 10 
and Hadamard waveform, the scaling laws also hold but with 
larger constants /j b- Therefore, depending on the number of 
receive antennas M r and the number of samples L in one 
pulse, mi and m 2 can be assumed to be bounded by small 
numbers. Based on [6 1, therefore, we conclude that the strong 
incoherence property is likely satisfied in our problem. 

It is also worth noting that under Gaussian orthogonal 
waveforms, 7772 is concentrated in a smaller range as com- 
pared with the range for Hadamard waveforms. This indi- 
cates that the waveform indeed plays a role for the use of ma- 
trix completion, and perhaps the waveform can be optimally 
designed to result in low probability of large values for m 2 . 

Since the matrix completion conditions appear to be sat- 
isfied in our case, we propose that each antenna obtains and 
forwards to the fusion center a small number of samples dur- 
ing each pulse. Note that along with the samples, the antenna 
needs to inform the fusion center on how the sample was ob- 
tained so that the fusion center can determine where to posi- 
tion the received samples in the receive data matrix Z q . In a 
practical setting, a pseudo-random sampling ADC at each re- 
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Fig. 3. Bounds of mi and 777,2 scales with reciprocal of 
and yL, respectively. 



Fig. 4. Performance comparisons: (a) Relative error of the 
recovered data matrix; (b) Probability of DOA estimation res- 
olution. 



ceive antenna could be used in the place of the Nyquist sam- 
pler, where the pseudo-random generator seeds would be dis- 
tributed to the receive antennas by the fusion center. Finally, 
the fusion center can recover the data matrix X 9 by applying 
matrix completion to the received data. 

Once the data matrix X g is recovered, it can go through a 
matched-filter bank to produce 



Y q = ^X q S H = B (9) SD,A T (9) + W„ 



(8) 



where W g is noise whose distribution is a function of the ad- 
ditive noise W 9 and the nuclear norm minimization problem 
in ©. Next, stacking Y q £ C M, xMt into a vector y q , and 
based on the vectors corresponding to Q pulses, the follow- 
ing matrix can be formed: Yr — [yi,...,yQ] G c MrMtX Q, 
Reshaping Y R into Y G C QM ' xMr , we have 



Y = FS[b (0 1 ) ) ... ) b(0*)]+W ) 



(9) 



where F = [d (t^) ® a (9 X ) , d (0 K ) <g) a (9 K )], d (0) = 
^ e j2n2-&T PRI ^ ^ e ji^{Q-\)T PRI ^ T ^ with g, denoting the 
Kronecker product. The sampled covariance matrix of the 
receive data signal can then be obtained as Ry = jj- YY H , 
based on which target estimation can be implemented using 
any array processing method such as MUSIC. 



5. NUMERICAL RESULTS 

In this section we present some simulation results on the per- 
formance of the proposed method. We use the simulation set- 
ting considered in the previous section, i.e., M t = 20, M r = 
40, Q = 5, L = 128. The signal-to-noise ratio (SNR) is set 
to 25 dB, while /3k are following complex Gaussian distribu- 
tion and kept unchanged for Q pulses. For matrix completion, 
the TFOCS software package [ 1 1 1 is used. 

First, we plot relative errors (averaged over 50 Monte 
Carlo runs) of the received data matrix X g for Hadamard 

and Gaussian orthogonal (G-Orth) waveforms. The relative 

|| z _x II 

error is defined as ,?„ M , where Z is the data matrix 
calculated without missing elements. Under each waveform, 



K = 2 point targets in the far field are randomly generated. 
The result is shown in Fig. |4](a) for q = 1. It can be seen 
from this figure that, as p increases, the relative recovery 
error of the data matrix under Gaussian orthogonal waveform 
reduces to the reciprocal of the SNR faster than that under 
Hadmard waveform. A plausible reason for this is that under 
the Gaussian orthogonal waveform, the maximum value of 
elements in the singular vectors of Z q is bounded by a smaller 
number with high probability, as compared with that under 
the Hadamard waveform (see Fig. |2). 

Next, the probabilities of DOA estimation resolution un- 
der the two orthogonal waveforms are plotted in Fig. |4](t>) 
for the following scenario. Two targets are randomly gen- 
erated among DOA range [—20°, 20°] with minimum DOA 
separations d B = [0.2°, 0.3°, 0.4°, 0.5°, 0.7° , 1°] and the cor- 
responding speeds are set to 150 and 400 m/s. The MUSIC 
algorithm is applied to obtain the target DOA information. If 



the DOA estimates 9; 



= 1, 



, 2 satisfy 

0. 1, we declare this as a success. The probability of DOA res- 
olution is then defined as the fraction of successful events in 
50 iterations. It can be seen from the figure that when p = 0.3, 
the Gaussian orthogonal waveform has a much better DOA 
estimation resolution compared with Hadamard waveform. 
As p increases to 0.5, the performance difference becomes 
small since the relative recovery errors under both waveforms 
are similar (see Fig. |4](a)). Figure |4] confirms that Gaussian 
orthogonal waveforms are better than Hadamard waveforms 
for matrix completion-based DOA estimation. 

6. CONCLUSIONS 

We have provided results suggesting that matrix completion 
can be used in MIMO radar to reduce the number of data 
needed to be communicated to the fusion center by each re- 
ceive antenna. Numerical results show that matrix comple- 
tion in conjunction with MUSIC can achieve accurate tar- 
get estimation with sub-Nyquist samples. Thus, the proposed 
method can result in significant savings in terms of data that 
need to be obtained at the receive antennas and subsequently 
transmitted to the fusion center. 
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